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We present an analysis of angular distributions and correlations of the X(3872) particle in the 
exclusive decay mode X(3872) — > J/%jj-K + n~ with J/ip — » We use 780 pb _1 of data from pp 

collisions at yfs — 1.96 TeV collected with the CDF II detector at the Fermilab Tevatron. We derive 
constraints on spin, parity, and charge conjugation parity of the X(3872) particle by comparing 
measured angular distributions of the decay products with predictions for different J PC hypotheses. 
The assignments J PC = 1 ++ and 2 _+ are the only ones consistent with the data. 

PACS numbers: 14.40.Gx, 12.39.Mk, 13.25.Gv 
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The recent discovery of the A (3872) particle 0, Q 
has revived general interest in charmonium spectroscopy. 
The exact nature of this particle is still unknown. At- 
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tempts to explain the A (3872) particle as a conventional 
bound quark-antiquark state have shortcomings, such as 
deviations from mass predictions or violation of isospin 
conservation [|. The close proximity of the A (3872) 
particle mass to the D°D*° mass threshold has raised 
the question whether the A(3872) particle is an exotic 
form of matter The determination of the quantum 
numbers spin J, parity P, and charge conjugation par- 
ity C is of vital importance for establishing the nature of 
the A (3872) particle. The evidence for the decay mode 
A(3872) — > J/ip 7 Q and the measurement of the di- 
pion mass distribution which is in agreement with 
the decay mode A (3872) — > J/ipp°, are consistent with 
a C-even assignment. Reference [6| observes an enhance- 
ment in the D°D°ir° mass spectrum and concludes that, 
if assigned to the A (3872) particle, low values for the 
spin quantum number are favored. Neglecting effects 
from model uncertainties in the dipion mass spectrum 
(see Q), preliminary results from Q favor J PC = 1 ++ . 
In this Letter we report the angular distributions in the 



decay A(3872) J/ip tt+tt", J/ip 



PC 



and com- 



pare them with predictions for different J states. The 
analysis is independent of any specific model of the in- 
ternal structure of the A (3872) particle. We consider 
all allowed states up to spin two and C-odd spin three 
states. 

We use a sample of pp collisions at a center-of-mass 
energy of ^fs — 1.96 TeV with an integrated luminos- 
ity of 780 pb _1 collected with the CDF II detector at 
the Fermilab Tevatron. The CDF II detector [ij consists 
of a magnetic spectrometer surrounded by electromag- 
netic and hadronic calorimeters and muon detectors. The 
tracking system is composed of a silicon micro-strip de- 
tector [9] surrounded by an open-cell drift chamber called 
the central outer tracker (COT) [lol l. We detect muons in 



planes of multiwire drift chambers in the pseudora- 
pidity range |r/| < 1.0. The J/ip — > fi + /i~ decays used in 
this analysis are recorded using a dimuon trigger, which 
requires two oppositely charged COT tracks matched to 
muon chamber track segments with an invariant mass 
from 2.7 to 4.0 GeV/c 2 . 

The basic event selection is described in although 
we do not cut on the dipion mass. Additional criteria 
are imposed on the number of candidates per event, the 
transverse momentum px of the A (3872) particle candi- 
date (> 6 GeV/c), the p T of the J/ip (> 4 GeV/c), and 
the kinetic energy released in the A (3872) particle decay, 
Q = m(J/ipTnr) - m(JHA - m{irir) (< 100 MeV/c 2 ), 
where m{J/ip) is from |l2j. The cuts are chosen to op- 
timize the significance S/\/ S + B of the observed signal, 
where S and B are the fitted number of signal and combi- 
natorial background events in a ±1.5<7 window centered 
on the A(3872) particle mass. The resulting distribution 
of the invariant J / 'tp ir + tt~ mass is shown in Fig. [TJ 

To simulate the decays of A (3872) particle states with 
specific J assumptions, we first generate phase-space 
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FIG. 1: The J/tp-K + ir~ mass spectrum after optimizing the 
selection cuts, fitted by a double Gaussian function for the 
ip(2S) (left), a Gaussian function for the X(3872) particle 
(right), and a second order polynomial for the combinatorial 
background. 



decays of X(3872) -> J/t/jn+w', J/ip -> n+ p,~ . De- 
tector effects are included using parameterized efficien- 
cies and acceptances. This sample is weighted according 
to each specific J PC hypothesis using the corresponding 
matrix element Mtot described below. 

The decay of the narrow X(3872) particle is mod- 
eled as the sequential two-body decay chain X(3872) — + 
J/ip (7r + 7r~), J/ip — > and the decay of the inter- 

mediate (7r + 7r~) state to 7r + 7r~. Assuming low relative 
angular momentum between the pions and conservation 
of C parity, the intermediate pion state can be in either a 
relative S'-wave ((7T7t)s) or a P-wave (p°) state. Mtot is 
formed by the product of a matrix element A4i for each 
decay and a term T(m(7T7r)), which describes the mass 
dependence of the intermediate (7r + 7r~) system. Because 
of the very narrow width of the intermediate J/ip, we can 
neglect the J/ip mass dependence. 

With fixed helicities the angular dependence of a two- 
body decay amplitude is given by the Wigner function 
D{\ Ai _ a . 2 [IE El, where Ji and \ are the spin and 
helicity of the decaying particle, and A^i and A^2 are the 
helicities of the child particles in the parent rest frame. 
The function is multiplied by two Clebsch-Gordan coef- 
ficients, coupling the spins of the child particles to their 
summed spin Si, and Si with their relative angular mo- 
mentum Li to Ji. 

In general, in the X(3872)— * J/ip (7r + 7r~) decay there 
is more than one combination to form J from L and 
S in a parity-conserving way. Of the independent am- 
plitudes corresponding to these combinations, only the 
ones with lowest L, assumed to be dominant, are taken 
into account. If more than one amplitude remains, mix- 
ing parameters are introduced to describe the physical 
state. Since the virtual photon in the J/ip — * de- 
cay can be treated as transverse, helicity combinations 



with 



A„- = are neglected. 
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The dependence of M.tot on the dipion mass has model 
ambiguities. Therefore, we do not use the information 
from m(7T7r) to distinguish between different J PC hy- 
potheses. The influence of the m(Tnr) model on the an- 
gular distributions via acceptance effects is very small. 
Nevertheless we choose for all J PC hypotheses the same 
model for the m(7T7r)-dependent terms, which agrees with 
the m(7T7r) spectrum measurement. In this way, no hy- 
pothesis is rejected due to a wrong m(7T7r) model. In 
detail, we fix T(m(irir)) to a relativistic Breit-Wigner for- 
mula with mass and width of a p° [H. Following @, we 
also fix the momentum dependence of the matrix element 
of the (7T+7T - ) — > 7r+7r~ decay to k* ■ fi (k*), where k* 
is the magnitude of the three-momentum of one of the 
pions in the (ir + ir~) rest frame and j\ (k*) is a Blatt- 
Weisskopf form factor 15[ to counter the divergence for 
rising k* . This form factor has the effective size r of the 
particle as a free parameter which we set to a common 
choice of r — 1 fm. 

A weight is formed from the square of the total matrix 
element M.tot by averaging over all initial state hclici- 
ties assuming unpolarized A (3872) particle production, 
incoherently summing over all final state helicities, and 
coherently summing over all intermediate state helicities. 

The decay is described by the decay angles 9x, Oj/tp, 
<Pj/i(i-> 4>Tnr, and A<I>, as defined in Fig. [2] For unpo- 
larized A (3872) particle production and because of ro- 
tational symmetry, the J PC of the A(3872) particle and 
the (ir + ir~) system affect the distribution of only four 
variables: m(7T7r), cos(6j/^), cos(0 7r7r ), and A<F 

The angular distributions are analyzed with a three- 
dimensional fit to take into account their correlations. 
From simulation studies, the optimal binning is deter- 
mined to be three bins in ||A<1> — w\ — ?|, and two bins 
in each of | cos(^j/^,)| and | cos^^)!, where absolute val- 
ues are used to exploit final state charge symmetry. The 
invariant J/-07r + 7r~ mass spectrum is fitted in each of 
the resulting 12 bins in a mass window of ±110 MeV/c 2 
around the A (3872) particle position using a binned max- 
imum likelihood fit, where the bin width is 2.5 MeV/c 2 . 
The distribution is described by a Gaussian function for 
the A (3872) particle and a second order polynomial for 
the background. The position and width of the Gaussian 
function describing the A (3872) particle are first deter- 
mined from a fit to the full invariant mass spectrum and 
are then fixed in the subsequent fits. We compare the 
fitted yield as a function of the angular variables with 
the predictions for different J PC assignments by forming 
a x 2 based on statistical uncertainties of the measure- 
ment. We determine the normalization of the simulated 
distributions from the measurement so that 1 1 degrees of 
freedom remain. 

The decay amplitude for the state with J PC = 1 h 
consists of three LiS-terms with the same L value; the 
J PC = 2 ^ state has an amplitude with two LS'-terms 
(see Table H]). None of the 1 h terms describes the data 



J PC 


decay 


LS 


X 2 (11 d.o.f 


) x 2 Prob. 


1++ 


J/ipp° 


01 


13.2 


0.28 


2"+ 


J/ipp° 


11,12 


13.6 


0.26 


\— 


J/-0(7r7r)s 


01 


35.1 


2.4xl0~ 4 


2+- 


J/i/>(7"0s 


11 


38.9 


5.5xl0~ 5 


1+- 


J/ip(ww)s 


11 


39.8 


3.8xl0~ 5 


2~ 


J/-0(7T7r) S 


21 


39.8 


3.8xl0~ 5 


3+- 


J/tJj(tttt) s 


31 


39.8 


3.8xl0~ 5 


3~ 


J/i/>(7"0s 


21 


41.0 


2.4xl0" 5 


2++ 


J/ipp° 


02 


43.0 


l.lxl0~ 5 


1 T 


ill 


10,11,12 


A C A 

45.4 


4.1x10 


o-+ 


J/ipp° 


11 


104 


3.5xl0~ 17 


o+- 


J/tP(tttt)s 


11 


129 


<lxl0- 20 


0++ 


j/ijp° 


00 


163 


<lxl0" 20 



TABLE I: Result of the X(3872) particle angular analysis. 
Listed are the state, the decay mode, the L and S quantum 
numbers of the J /ip-(ir + -K~) system, the \ 2 with 11 degrees 
of freedom and the x 2 probability. 



alone, so we fit for a mixed state by minimizing the x 2 - 
For the 2 h state, the amplitude for S = 1 is sufficient 
to describe the data. 

Table U shows the x 2 f° r each J PC assignment. We 
find that only the assignments J PC — 1 ++ and 2 ^ are 
able to describe the data. All other states are rejected 
by more than 3 standard deviations (x 2 prob. < 2.7 x 
10~ 3 ). Figure[3]shows the measurement and the expected 
distribution for four of the assignments. 

An important cross-check of the analysis is to verify 
whether the correct result is obtained for the ip(2S), 
with known quantum numbers J PC = 1 , which de- 
cays into the same exclusive final state as the A (3872) 
particle. For the 1 assignment, the fit probability is 
1.5%. Using the ip(2S) model of Novikov and Shifman 
[l6| . which includes a small D-w&ve admixture in the 
description of the (tt + it~) system, the fit probability is 
17.9%. The sensitivity to such a small admixture is only 
present in the high statistics ip(2S) sample. The next 
best model J PC = 2 ++ has a fit probability of 0.58%, and 
all other hypotheses that were tested yielded fit proba- 
bilities smaller than 2 x 10~ 6 . 

We vary several inputs to the fitting procedure and 
the model of the A (3872) particle to investigate the sta- 
bility of the x 2 ■ Figure 0] shows the resulting x 2 values 
for the different J PC hypotheses for the variations in- 
vestigated. The default analysis is shown as variation 
(1). The following effects arc considered: (2)/(3) de- 
crease/increase the fit window by 20 MeV/c 2 , (4)/(5) 
decrease/increase the bin width to 2.0/2.86 MeV/c 2 , 
(6)/(7) vary fixed A(3872) particle mass by ±lcr, (8)/(9) 
vary fixed X(3872) particle width by ±lcr. 
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cos(e. j| > 0.6 



coste^i 

>0.5 




0.63 1.15 %I2 0.63 1.15 7t/2 

0.63 1.15 ti/2 0.63 1.15 jc/2 

||A<J>-7t| - tt/2| 

FIG. 3: Measured 3D angular distribution with acceptance 
corrected predictions for J PC = ++ (solid line), 1 ++ (dotted 
line), 2~ + (dashed line), and 1 (dash-dotted line). The 
plot is divided into 2x2 regions, corresponding to intervals of 
|cos(^j/^)| and | cos(^ 7r7r ) | ■ Each region shows the distribu- 
tion of ||A$ — 7r| — -|| in 3 bins. The bin contents have been 
scaled to the same bin volume. 
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FIG. 4: Total x 2 f° r different analysis variations on the y- 
axis, explained in the text. Vertical bars are added for visual 
guidance. The \ 2 values of the spin states are all above 100. 
The 2 and 3" 1 states have the same angular distribution 
as the l +_ state. 



To evaluate the contribution to the systematic uncer- 
tainty from our choice of the m(inr) spectrum, the fol- 
lowing variations are considered: (10) fix form-factor r 
to 0.001 fm, (11) fix form-factor r to 100.0 fm, (12) use 
simple phase space for m(7T7r). 

Finally, systematic uncertainty due to details concern- 
ing the simulation has been considered by varying distri- 
butions for (13) p T and (14) r\ of the A(3872) particle, 
switching off (15) a. pt dependent efficiency correction 
for the pions, (16) a <f> dependent correction of the COT, 
and (17) an effective rj correction used to model the po- 
sition of the generated primary vertex. All variations 
are consistent with 1 ++ and 2 h being the only likely 
assignments. 

A conventional explanation for the A(3872) resonance 
is a charmonium (cc) state. In this picture, the state 



with J PC = 1 ++ could be identified with the x'd an( i the 
assignment J PC = 2 ^ with the r/ C 2. An exotic interpre- 
tation is that the A (3872) particle is a molecular state 
or that a significant four-quark interaction contributes 
to the wave- function [l7|. The result of this analysis is 
compatible with the models of a molecular state devel- 
oped by Tornqvist [l8| and Swanson (lij , who predict the 
quantum numbers J PC = 1 ++ for a bound D D* state. 

In summary, a spin-parity analysis of the A (3872) par- 
ticle in the final state 7r + 7r~ has been performed. 
The method of helicity amplitudes has been used to ana- 
lyze A(3872) -> J/xp(irir) s and A(3872) -► J/ifjp tran- 
sitions. Using a x 2 approach to compare expected angu- 
lar distributions with measured distributions, it is found 
that only the C-even assignments J PC = 1 ++ and 2 
both decaying via J/ifjp°, describe the data. All other 
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states are excluded at 99.7% confidence level. 
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